Periodically driven Quantum Ratchets: Symmetries and Resonances 
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We study the quantum version of a tilting and flashing Hamiltonian ratchets, consisting of a 
periodic potential and a time-periodic driving field. The system dynamics is governed by a Floquet 
evolution matrix bearing the symmetry of the corresponding Hamiltonian. The dc-current appears 
due to the desymmetrization of Floquet eigenstates, which become transporting when all the rel- 
evant symmetries are violated. Those eigenstates which mostly contribute to a directed transport 
reside in phase space regions corresponding to classical resonances. Quantum dynamics leads to 
the dependence of the average velocity on the initial phase of the ac-field. A resonant enhancement 
(or suppression) of the dc-current, due to avoided crossings between different Floquet states takes 
place upon tuning some control parameters. Our studies are predominantly aimed at experimental 
realizations of ac-driven quantum ratchets with cold atoms. 

PACS numbers: 05.45.Mt, 05.60.-k, 32.80.Pj 
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I. INTRODUCTION 

Ratchets are viewed as realizations of systems which 
produce a directed current from a fluctuating envi- 
ronment in the absence of gradients and net forces 
[H, 0, i, i, ||. Initially proposed as an abstract phys- 
ical model for the understanding of a micromolccular 
machinery ratchet systems have found diverse ap- 
plications in many areas 0, 0, [E], from mechanical de- 
vices jgj up to the quantum systems @, H, H, [HJ [HI- 
Among the areas with growing interest is the study of 
ratchet dynamics for cold atoms [l2j . Hamiltonian ratch- 
ets [III, 0, fl5l Jl6l . H3] with the corresponding symmetry 
predictions [13| have recently been successfully realized 
with cold Rubidium and Cesium atoms in optical lattices 
with a two-harmonics driving and a tunable weak dissi- 
pation [18j . In these experiments, the mechanism of the 
Sisyphus cooling [l9| has been used in order to furnish 
initial conditions in form of an optical lattice: an ensem- 
ble of atoms localized in the wells of a periodic potential 
created by laser beams. In the momentum space this 
corresponds to a narrow distribution near the momen- 
tum p = 0. This is essential for the observation of the 
rectification effect, since the asymptotic current tends to 
zero for broad initial distributions in momentum space. 

To describe the ratchet dynamics of the thermal cloud 
of cold atoms in optical lattices [l8[ , a classical model has 
been used [HI, [lj, [HI]. However, as the cloud of atoms 
in a far detuned optical potential [2(| is cooled further 
down, one can obtain a system where quantum effects 
become relevant. Many studies of quantum ratchets are 
based on the kicked rotor model @, [l6|, [2l|, \2^, which 
is readily treated numerically. In reality kicks may heat 
the system and degrade quantum effects. It is perhaps 
for that reason, that other successful experiments [Hj] use 
a time continuous drive, which is also the choice in the 
present work. 

The paper is organized as follows. In Section II we 
introduce the model for the quantum tilting ratchet and 



discuss its properties. The dynamics of the system is 
studied within the Floquet theory. We use the Husimi 
representation in order to link eigenstates to manifolds in 
the classical phase space. In section III we study the rel- 
evant symmetries, whose violation leads to the directed 
current, and the desymmetrization mechanisms. Section 
IV is devoted to the dynamics of current rectification. 
We begin by discussing the dependence of the current on 
the initial condition. Then we investigate the appear- 
ance and controlling of quantum resonances that lead to 
a significant enhancement of directed transport. In Sec- 
tion V we introduce the realization of a quantum flashing 
ratchet system. Finally, in Section VI we summarize our 
results. Some of our important technical details are de- 
ferred to the Appendices A, B and C. 



II. MODEL SET-UP 

Let us start with a cloud of atoms, placed into a peri- 
odic potential (formed by an optical standing wave) and 
exposed to an external ac field. Assuming a low density, 
we neglect the interactions among the atoms. Hence, our 
problem can be described by the Schrodinger equation 



m-m)) = m,to)m)), (i) 

where the Hamiltonian H is of the form [l8j] 

H(x,p,t) = ^- + (l + cos(x))-xE(t-t ). (2) 

Here E(t) is an external periodic field of zero mean, E(t — 
t a + T) = E(t - t ), J T E(t)dt = 0, and t is the initial 
time. 

The Hamiltonian ([2]) is periodic in time with pe- 
riod T. The solutions \ip a (t)) = U(t,t )\t(j a (0)) of the 
Schrodinger equation (1) can be characterized by the 
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eigenfunctions of U(to +T, to) which satisfy the Floquct 
theorem: \ij) a {t)) = e^*^)), |0 a (t + T)) = |0 a (t)). 
The quasienergies £" Q (— 7r < £? Q < 7r) and the Floquct 
eigenstates can be obtained as solutions of the eigenvalue 
problem of the Floquet operator 

U{T,t )\Mto)) =e- iE -\Mto))- (3) 

The Floquet eigenstates provide a complete orthonor- 
mal basis and the stroboscopic quantum state can be 
expressed as [23| 

\iP(mT,t )) =J2 C «(to)e- mE «\Mto)), (4) 

a 

where the coefficients {C a } depend on to. For later 
convenience the integer a = 0,1,2.., sorts the states 
\4> a ) such, that the mean kinetic energy (p 2 ) a = 
1/T J Q (4> a \p 2 \<Pa)dto monotonically increases. 

By using the gauge transformation, \ip) — > 
cxp(— J-x J* E{t')dt')\ip) [24|, we transform the original 
Hamiltonian {(2j) into a spatially periodic one (see 
Appendix A). Then the solution of the time-dependent 
Schrodinger equation for the new Hamiltonian may be 
written as 

\lp(t)) = e -kfo{klP~Mt',t„)]* + (l+cosx)}dt' |^ ^ ^ 

with the vector potential A(t) = — f Q E(t')dt'. Due to 
discrete translational invariance and Bloch's theorem all 
Floquet states are characterized by a quasimomentum k 
with \<f> a (x + 2tt)) = e lhK \(/) a (x)). 

Here we choose n = which corresponds to initial 
states where atoms equally populate all (or many) wells 
of the spatial potential. The wave function is expanded 
in the plane wave eigenbasis of the momentum operator 



p, \n) 




x 

\m) = E c «wi™>- ( 6 ) 

n——N 

A detailed discussion of the two numerical procedures 
used to integrate Eq.([T]) is presented in the Appendix IBl 

To examine the morphology of the quasienergy states 
we use the Husimi representation [25| 

p(( x ) 7 {p)) = -L m <S, {xh{p) )f, (7) 

where 

$ (->.(p)W = (2^2)1/4 eX P ( 4^2— +*<P> S j 

(8) 

with a — (h/2) 1 / 2 . Here (x) and (p) stand for the average 
of the position and momentum respectively. The Husimi 
representation provides an insightful, coarse grained de- 
scription in the phase-space [251 ]. 



III. RELEVANT SYMMETRIES AND THEIR 
VIOLATIONS 

A. Classical limit 

Let us briefly outline the classical case. In this limit the 
system is generically characterized by a nonintegrable dy- 
namics with a mixed phase space containing both chaotic 
and regular areas in the three-dimensional phase space 
(x,p, t) [13] • Due to time and space periodicity of the 
classical equation of motion x = — we can map 
the original three-dimensional phase space onto a two- 
dimensional cylinder, T 2 = (x modlip), by using the stro- 
boscopic Poincare section after each period T = 27r/w, 
cf. Fig. 1(a). This provides a helpful visualization of the 
mixed phase space structure. A stochastic layer, located 
near the line p — 0, originates from the destroyed separa- 
trix of the undriven system, E{t) = 0. The chaotic layer 
is confined by transporting KAM — tori, which originate 
from perturbed trajectories of particles with large kinetic 
energies, (p 2 )- The stochastic layer is not uniform and 
contains different regular invariant manifolds - regular 
islands, unstable periodic orbits and cantori [27l |. 

There are two symmetries which need to be broken 
to fulfill the necessary conditions for the appearance of 
a dc-current J = lim t ^ 00 | = (x) [HI]. If E(t) is shift 
symmetric, i. e. E(t) = —E(t + T/2), the symmetry 

S a :(x,p,t)^(-x,-p,t + T/2), (9) 

is realized. Furthermore, if E(t) is symmetric, i. e. 
E(t + t s ) = E{—t + t s ) at some appropriate point t s , 
the symmetry 

S b : (x,p,t) -> (x,-p,-t + 2t s ), (10) 

is holds. Any trajectory, when transformed using one of 
the abovementioned symmetry operations, yields again 
a trajectory of the system but with opposite velocity. 
Assuming that ergodicity holds in the stochastic layer 
(which means that all its average characteristics are in- 
dependent on the initial conditions, and that symmetry- 
related trajectories have the same statistical weight, we 
conclude that the asymptotic velocity within the chaotic 
layer is zero. So, whenever S a and/or Sb are realized, 
directed transport is forbidden inside the chaotic layer 
0. 



The two frequency driving, 

E(t,t ) = E 1 cos[u(t~t Q )}+E 2 cos[2uj(t~to) + 9}, (11) 

ensures that for E\,E 2 ^ S a is always violated. In 
addition Sb is violated for 9^0, ±tt. The appearance of 
a nonzero dc-current J c h — lim^oo 1/t J t p(t')dt' in this 
case is due to a desymmetrization of the chaotic layer 
structure (FigJTJi) . It induces a desymmetrization of the 
events of directed motion to the right and left [28( . Due 



FIG. 1: (a) Poincare section for the classical limit; (b-d) 
Husimi representations for different Floquet eigenstates for 
the Hamiltonian ([2]) with h = 0.2 (momentum is in units of 
the recoil momentum, p r = hkn, with kh = 1). The parame- 
ters are E\ — E2 = 2, lo = 2, 8 — —tt/2 and to — 0. 

to ergodicity inside the layer, the asymptotic current is 
independent of the initial time to, for initial conditions 
located inside the chaotic layer. With the specific choice 
of the driving (|TT|) it follows J c h(6) — —Jch(-O) and 
J c h(9) = —Jch{8 + 1") (l3| . From perturbation theory 
it follows J c h ~ EfE 2 sin 

El. An efficient sum 
rule allows to compute the average current J c h by proper 
integration over the chaotic layer [l6| . 

B. Quantum case 

If the Hamiltonian is invariant under the shift sym- 
metry S a ©, then the Floquet operator possesses the 
property, see Eq.(179) in Ref.[H|: 

U(T, t ) = £7* (T/2, t f U (T/2, t ) . (12) 

Here £7* performs a transposition along the codiagonal 
of U. With (Tnj) 5 a is always violated. 

Likewise, one can show that, if the Hamiltonian is in- 
variant under the time reversal symmetry St, (|10p . then 
the Floquet matrix has the property [29(] 

U(T,t ) = U(T,t f (13) 

That symmetry will be recovered for = 0, tt. More 
details about the derivation of the previous properties of 
the operator evolution are explained in the Appendix [C] 

We introduce the quantum asymptotic current as the 
mean momentum expectation value 

■/(to) - (P)t = Urn l/t f dt' (iP(t',t )\p\iP(t',to)), (14) 

where p = (ip\p\tp) is the mean instantaneous momen- 
tum. Expanding the wave function in the Floquet basis 
the current becomes 

Ak) -^(p)a|C Q (i )| 2 , (15) 
a 
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FIG. 2: A part of the quasienergy spectrum as a function 
of the parameter 9. The symbols indicate the corresponding 
Floquet states shown in FigfIJb-d). The empty red circle 
indicates the existence of an avoided crossing between two 
eigenstates. The empty square enclose the region of spectrum 
that appears enlarged in Fig(3] 



where (p) a is the mean momentum of the Floquet state 
\4> a ) (see Appendix |D] for more details). Thus, we have 
to study the properties of Floquet states. 

In general, the Floquet bands of the system (JlJ^J) de- 
pend on the control parameters like amplitudes, frequen- 
cies and phase shifts of the components of the force E(i). 
We vary the parameter while keeping the other fixed. 
In FigJ5] we present the quasienergy bands as a func- 
tion of 0. There are two remarkable features. First, the 
spectrum exhibits two symmetries, E a (0) = E a {—9) and 
E a (0) = E a (0 + ir), which are consequences of the choice 
(TrTj) . Second, while some bands show strong dispersion 
upon the variation of 0, reaching the maximal dispersion 
at = ±7r/2, others have a rather flat dependence. In 
FigIHb-d) we present Husimi functions [25| for several 
Floquet states depicted by symbols in FigO The Planck 
constant h — 0.2 is in a range, where it is possible to es- 
tablish a correspondence between different Floquet states 
and the invariant manifolds of the mixed phase space for 
the classical limit. The states (b-d) are located in various 
regular phase space regions. 

On the other hand, from Figf2]we observe repulsion be- 
tween some bands whereas other apparently cross each 
others. These crossings are in fact avoided crossings 
which lead to interesting effects. 

For the symmetric case = 0, ±7r the Floquet ma- 
trix has an irreducible representation using even and odd 
basis states \n) Sia = (\n) ± | — n))/ y/2 — S^q. All the 
Floquet states appear as quasidegenerated doublets and 
(p) a = for all a (see FigsHHSJ). We especially note that 
this is true for states with arbitrarily large kinetic en- 
ergy, for which the corresponding quasienergies become 
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almost pairwise degenerated. Consequently, following 
the Eq. (fT5|) , J = in this case. 

For 9 7^ 0, ±7r the Floquet states become asymmet- 
ric and the quasidegeneracies are removed (Fig[3]). Flo- 
quet states acquire a non-zero mean momentum, thus 
becoming transporting. On the other hand, while in the 
symmetric case it is possible to have coherent tunnel- 
ing oscillations between disconnected regular islands, in 
the non-symmetric case a possible dynamical tunneling 
is suppressed [3(|. Thus, a wave packet initially located 
on a Floquet state with an asymmetric distribution of 
momentum (similar to the state shown in FigOJ)) will 
undergo a permanent directional transport. 

This phenomenon of desymmetrization results in a non 
trivial dependence of the momentum upon 9. We can 
gain some understanding of this effect by modeling the 
evolution of two eigenstates, which form a parity-related 
pair at the point 9 = 0, by neglecting their coupling with 
the rest of the states. 

To this end, we take two eigenstates, $ a and which 
can be expanded in the plane wave basis ±|n). The cor- 
responding propagator for such a reduced system is 



U ab 



A 



A 
1 + e 



(16) 



where the nondiagonal term A refers to the interac- 
tion between the states. In general, this quantity de- 
pends upon the momentum and driving field's charac- 
teristics, i.e. A = A(p, Ei, E 2 , 9). On the other hand, 
the parameter e breaks the parity symmetry, and yields 
e ~ nE\E 2 sin0 at a first-order of a perturbation theory. 
It is known that, in the absence of driving, E(t) = 0, 
the asymptotic dependence of the splitting value on n is 
A ~ nT n , i.e. a superexponential decay [3l|. 
The eigenstates in a plane- wave basis are 



4\, 



1 



i 

W) 



L • |n) + (7 + Vl + 7 2 ) • I - n>] , (17) 
(7 + Vl + 7 2 ) • \n) - 1 • I - n>] , (18) 



where 7 = e/A and L{i) = \jl + (7 + y/l + 7 2 ) _ . The 

symmetric case, 9 = 0, corresponds to the limit 7^0. 
The eigenstate's momenta are 



Ps 



~in, 



Pa « in, 



At the opposite limit, 7 — > 00, we have 



Ps 



-n 1 



1 

472 



p a K, n 1 + 



47 5 



(19) 



(20) 



For large n, A is small, but for e = (symmetric case) 
it follows 7 = 0, and the eigenstates carry zero momen- 
tum. Deviating from 9 = 0, the moments quickly reach 
values ±n. That happens for E\E 2 sin# 3> n - ™ -1 , which 
will hold for any nonzero E\, E 2 , 9, provided n is large 
enough. 




FIG. 3: Right panel: Quasienergy spectrum in the vicinity of 
two quasidenerated Floquet states. The symmetry is restored 
at the points 9 = 0, ±7r where the splitting is rather small. 
Left panels: Husimi function for different 8 a) 9 = —0.05; b) 
6> = 0. 



Leaving the region corresponding to the chaotic layer, 
the splitting between the doublets drops quite suddenly 
to zero. This explains the momentum asymmetry in Fig [3] 
for values of 9 slightly deviated from 9 = 0. 



AVERAGE CURRENT AND RESONANCES 
IN CURRENT RECTIFICATION 



IV. 



The asymptotic current fTJJ depends on the initial 
conditions, and, following Eq.fTS]), only those Floquet 
states which overlap with the initial wave function, ip(to), 
contribute to its subsequent time evolution. From now 
on, we restrict our analysis to the initial condition in 
the form of a plane wave with wavevector n = 0, i.e. 
ip(to) = |0) = This initial condition spreads over 

Floquet states with low kinetic energies, namely over 
those eigenstates for which the Husimi representations 
lie in the chaotic layer's region (Fig|4]). 



In the previous section we have discussed the mech- 
anism of the dessymetrization when we tuned the sys- 
tem away from the symmetry point 9 = 0. In or- 
der to estimate a dc-current value we need a quali- 
tative information about mean momentum values ac- 
quired by Floquet states. As a suitable quantity to 
study this issue, we use the cumulative average momen- 
tum, P a = J2 s <a(p) a - FigH] depicts the cumulative av- 
erage momentum, obtained from the recursive relation 

Pa+2 =Pa + (p) a +l + (p) a +2, with P = (p) [26]. The 
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FIG. 4: Top: The line corresponds to the cumulative momen- 
tum, P a (see text). All the momenta are scaled in units of the 
recoil momentum p r . Inset: The bars depict the weights of 
the respective Floquet states that overlap with the initial state 
jO). The bars have been multiplied by a factor of 10. Bottom: 
The plot of the average kinetic energy, complementary to the 
Poincare section on Fig. 1(a). Starting from an initial point 
(x(to),p(to), we calculate the kinetic energy averaged over one 
period, E kin (mT) = ^ f*°^^ r _ 1 ^ T p 2 (t)dt, and consequently 
plot points (x(mT) , Eki n (mT)) . The parameters are the same 
as in Fig[T] 



asymmetry stems mainly from Floquet states located 
in the chaotic layer region of the classical phase space. 
The dependence has several peaks. The Floquet states 
with higher values of the mean momentum, which yields 
strong dc-currents, are located at kinetic energy's regions 
that correspond to classical transporting resonance is- 
lands (see bottom part of the Fig0]). A more detailed 
analysis by using the Husimi distribution confirms this 
finding. There are infinitely many high-order resonance 
islands in the phase space of the Hamiltonian system (2) 
[27l | . whose sizes shrink quickly as the resonance's order 
increases with (p 2 ) — > 00. The desymmetrization effect is 
observable only when a resonance island can host more 
than one Floquet state. With increasing {p 2 ) a , P a ap- 
proaches zero in full accordance with the fact that total 
current over the whole momentum space should be zero 

The asymptotic current for a fixed initial condition 
\ip(t(j)) depends in general on the initial time to- Note 
that this is possible also in the classical case, since the 
initial distribution may overlap with different regular 
transporting manifolds ,16]. However, if we start with 
a cloud of particles located exactly inside the chaotic 
layer, the classical asymptotic current will be indepen- 
dent of to for any choice of the distribution function over 
the chaotic manifold. This is not true for the quan- 
tum case where the current may even change its sign 
with the variation of t . It is a consequence of the 
linear character of the Schrodinger equation 32]. We 
will first discuss the results obtained after averaging over 
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FIG. 5: The average current J (in units of the recoil momen- 
tum) vs 9. Blue and red delta points are obtained with higher 
resolution in particular ^-regions. Insets: (top right): En- 
largement of the region around 9 = —2.06. (bottom): Husimi 
function for the eigenstates in the avoided crossing indicated 
in FigH 



the initial time t . Then we can assign a unique cur- 
rent value, J = l/T J(t )dto, for fixed parameters 
of the ac-field, E\, E2, and 9. Figj5] shows the depen- 
dence of the average current on the asymmetry param- 
eter 9 for the initial condition = |0). The aver- 
age current J shows the expected symmetry properties 
J{9) = -J(8 + 7r) = -J(-0). These symmetries for the 
current imply that the results obtained for the interval 
(— 7T, —tt/2) holds for other intervals as well. 



The typical dependence of the average current J on 
9 is shown on Fig(5]The curve presents a smooth profile 
with several peaks. 

By comparing Figj5] with the quasienergy spectrum 
FigEl one can associate peaks with single sharp avoided 
crossings (resonances) between two Floquet eigenstates. 

The Husimi function gives us additional insights for 
the appearance of these peaks (see insets in Fig. 5). It 
shows two states that are located in the chaotic layer, 
one of them near p = 0, while the other is off the line 
p = with a strong asymmetry in momentum, corre- 
sponding to a regular transporting manifold. Off reso- 
nance the initial zero-plane wave state mainly overlaps 
with the first eigenstate, which yields some nonzero cur- 
rent due to desymmetrization. At resonance this eigen- 
state mixes with the "transporting" one, resulting in a 
strong enhancement of the current which is reflected in 
the presence of peaks. 

From an experimental point of view a too narrow res- 
onance may become undetectable due to resolution lim- 
itations. We thus studied how to vary the width of the 
resonance without much affecting its amplitude. It turns 
out to be possible by tuning another control parameter, 
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FIG. 6: (a) The average current J (in units of the recoil 
momentum) vs 6 for different amplitude values of the second 
harmonic, Ei: 0.95 (pointed line), 1 (dashed line) and 1.2 
(solid line). Insets: relevant details of the quasienergy spec- 
trum versus 9 in the resonance region for E% = 1 (top right) 
and E2 = 1.2 (bottom left). The parameters are E\ = 3.26 
and uj — 3. (b) Husimi function for the upper eigenstate that 
appears in (a) (top right inset) with 8 = —n/2. 



e.g. the amplitude E 2 - We decrease this field amplitude 
in order to disentangle the two Floquet states and remove 
the avoided crossing. That will happen for some value of 
E 2 at 9 = ±7r/2. The details of the quasienergy spec- 
trum around that critical point are shown in the insets 
in Figj6£i. The two quasienergy spectra disentangle for 
E2 = 1 but stay close over a sufficient broader range of 9 
values. Thus the resonances become broader, as seen in 
FigJSJi. Further decrease of E2 to a value of 0.95 leads to 
a strong separation of the two spectra, and consequently 
to a fast decay of the amplitude of the resonance. 

As we already know, in resonance Floquet states mix. 
Here the new eigenstate contain contributions both from 
the original chaotic state as well as from the regular 
transporting island state (see FigJS]). The Husimi dis- 
tribution of this mixed state shows strong asymmetry as 
expected. The states from the chaotic sea and the regular 
transporting island are connected by a narrow isthmus. 
Therefore, it is possible to leak from the nonstransport- 
ing to the transporting island, which, in principle, should 
reduce significantly the observation time for the detec- 
tion of resonances. It was shown in [28[, for the classical 
limit, that for a maximum desymmetrization of the phase 
space, particles may stick to a regular island leading to 
a ballistic flight. Hence, based on the classical analog, 
we can say that our hybrid state serves as a quantum 
ballistic channel. 
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FIG. 7: Current dependence on the initial time to and 9. The 
parameters are the same as in Fig(6] 

On the other hand, as already mentioned, the asymp- 
totic current depends on the initial time to- The ob- 
served resonance structures, due to resonant interaction 
between Floquet states (avoided crossings of quasiener- 
gies), are independent of the initial time to. In Fig[7]we 
plot the nonaveraged current as a function of both 9 and 
to- While the smooth background is barely resolvable 
with the naked eye, the resonances are clearly seen, and 
their position is not depending on to, while their ampli- 
tude does. That implies that one can further maximize 
the resonant current by choosing properly the initial val- 
ues for to- 

So far, we have considered only the frequencies, to = 
2, 3 (Figs. 1-7). To gain a better understanding of the res- 
onances, we compute the dependence of the current as a 
function of the frequency of the driving force (see Fig [8]). 
This figure shows the appearance of peaks, whose ampli- 
tudes increase as the frequency decreases. However, for 
larger frequencies, peaks show up for specific frequencies, 
a result which depends on the drive amplitudes. Inter- 
estingly, one can notice the inversion of current. 

In order to increase the current, we reduce the fre- 
quency of the ac- force E(t). The asymptotic current for 
uj = 1 is shown in the inset of FigJS] We again observe 
a clear broad resonance, but the maximum current value 
increases by an order of magnitude up to 0.5 in units of 
recoil momenta. 

Real experiments are always limited by a finite obser- 
vation time. Therefore, the rate of convergence to the 
asymptotic current value, Eq.(12), becomes a crucial is- 
sue. To give a more realistic approach of our above results 
to experiments, we perform numerical computations us- 
ing the running average momentum (current) definition, 

1 r* 

P=- Pdt. (21) 
1 Jo 

The time evolution for the instantaneous momentum 
P(t) is shown in FigJSJ Top of FigJS] shows that off reso- 
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FIG. 8: The average of current J versus frequency uj of the 
driving force for 9 — —n/2. Inset: The average of current J 
versus 9 for w = 1. The plot of the current is depicted for 
the interval (— 7r/2, 0), since the current is antisymmetric with 
respect to a reflection at the origin. The other parameters are 
E x = 3, E 2 = 1.5. 



nance the momentum rapidly evolves towards an asymp- 
totic value, while in resonance the momentum performs 
long term oscillations which slowly approach the asymp- 
totic value as the time increases. To speed up the conver- 
gence of the momentum towards the asymptotic value, we 
reduce the frequency. Figj9] shows large currents close to 
the recoil momentum with relaxation times of the order 
of 100 to 200 periods for u> — 1. The above results were 
obtained for h — 0.2. We have repeated the calculations 
using h ~ 1 and observed similar results. 
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FIG. 9: Momentum evolution of state |0) in units of periods. 
Top panel: On resonance cases 9 = —n/2 for different to. Off 
resonance case 9 = —2.2. The parameters are u> = 3, E\ = 
3.26, E2 = 1. Bottom panel: Case ui — 1 with parameters 
£1 = 3, E 2 = 1.5. 



V. QUANTUM HAMILTONIAN RATCHETS 
WITH FLASHING POTENTIALS 



Recent experiments have shown the possibility to 
achieve an optical lattice with variable asymmetry [33j . 
The Hamiltonian in this case is given by 



S 1 :(x,p,t)->{-x + 2x.,-p,t) (24) 
if U(x + x s ) = U {—x + x s ) for some x s , 



H(x,p,t) 



+ U(x)E(t-t ), 



where U(x + L) = U(x) and E(t + T)= E(t). 
Here, we take, as in [33j . the potential 



U(x) = K [cos{x) + s cos(2x + 9 P )} , 



(22) 



(23) 



where is the parameter which makes the potential 
asymmetric for values different from 0, ±7r. The relevant 
symmetries for the classical limit of the Hamiltonian 
are 



S 2 : (x,p,t) -> (x,-p,-t + 2t s ) (25) 
if E(t + t a ) = E(-t + t s ) for some t s , 



S 3 : {x,p,t) -> (-x + 2x s ,-p,t + T/2) (26) 
if E(t) = -E(t + T/2) and U(x + x s ) = -U(-x + x s ), 



S 4 : (x,p,t) -» {x + L/2,-p,-t+2t a ) (27) 
if E(t + t s ) = -E(-t + t s ) and U(x) = -U{x + L/2). 
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FIG. 10: The average current J vs 6. Inset: Asymptotic 
current as a function of the initial condition to at 8 = —it/2. 
The parameters are E\ = 2, U2 = 1.5, ft = 1, w = 1, if = 1.5, 
s = 0.25, pot = vr/2. 



It is easy to see that a two-harmonic potential Eq. 
is insufficient to break the time-reversal symmetry S2 
1(250. Thus an asymmetric ac- modulated function E(t) 
is needed to desymmetrize the system. We use the 
two-harmonic function (|11|) , which ensures that for 
p ^ 0, ±7r and 6* 7^ 0, ±7r all the relevant symmetries, 
Eqs. (|24l25l26l27p . are violated. As stated before the 
quantum system inherits the same symmetries as the 
classical limit. Therefore a current should appear as the 
symmetries are broken. 

FigHn] shows the average current J versus 9 for the 
asymmetric potential, 9 p = — tt/2, at strong quantum 
limit H = 1. Since the current dependence possesses the 
symmetry J(9) — — J(9 + n) = —J(—9), it is sufficient to 
plot the current for the interval 9 e (— 7T, 0). As for the 
case of rocking ratchets, we found evidence for quantum 
resonances and dependence of the asymptotic current on 
the initial time to (see inset in FigfTT 



VI. CONCLUDING REMARKS 

In summary, we have investigated the mechanisms 
of the emergence of an average current in ac-driven 
quantum systems. The key mechanism of the current 
rectification is the desymmetrization of Floquet states. 
A peculiarity of the quantum ratchet is its strong 
dependence of the current on the chosen initial time to 
of the applied field. We identified quantum resonances 
induced by avoided crossings between Floquet states 
which enhance the current drastically. Optimizing the 
drive frequency, amplitude and initial phase, resonant 
currents easily reach the recoil momentum value and 
should be experimentall y o bservable using driven cold 
atoms in optical lattices [l8l |33| . 
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APPENDIX A: GAUGE TRANSFORMATION: 
SOME RELEVANT QUANTITIES 



We introduce the transformation 

i 



\tp) = exp 



-xA(t) 



with A(t) = - f*E(t')dt'. 

Substituting into Eq.([T]) we can rewrite the 
Schrodinger equation as 



(Al) 



where H(t) = i 



AU) 

p--r 



+ [1 + cos(x)]. 

To compute the momentum we start from the defini- 
tion 

d 



p(t) = &\-m—\i)) = 

f°° ~ d ~ 
-ih tp* —ipdx = p(t) + A(t) 
J -00 ox 

where (p) = (p). 

The kinetic energy T is given by 



(A2) 



T(t) = ( i> 



exp 



-xA{t) 



if, 



dx 2 



dx 2 



exp 



4> 



--xA(t) 



After straightforward calculations 

T(t) = f[t) - 2A(t)f(t) + A(t) 2 . 
The average kinetic energy reads 

(T) = {f)-2{A(t)f(t)) + {A(t) 2 ). 
Recalling the Husimi distribution, 

p((x),(p)) = ^\{il>\9 {at)M )\ 2 , 
we can recast it in the new frame of variables as, 



V- ) ■ (A3) 

(A4) 
(A5) 



(A6) 



P((s}i <P» 



1 

2V 



e" 4 1 Mt) ijj* (x, t ) e l {p) x+ Li ^ )1 dx 



1 

2V 



rixM {p) * ]x+ -^^dx 



(A7) 
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Finally, after reversing the operations, we find 



we obtain for 4> n (t) the solution 



... ;( X ) ; ( P ) + MJ 



(A8) 



APPENDIX B: NUMERICAL INTEGRATION 
SCHEMES 

In order to find the numerical solution for the time- 
dependent Schrodinger equation we have used two differ- 
ent methods of integration. The first method has been 
described in (34j . It starts by using the expansion Eq. © 
and approximating 



A(t) = ^A(t)Y / S(t-kAt), 



(Bl) 



for the time dependence of the vector potential A(t) = 
- f* o E{t')dt', with At = 2n/(ujN t ), where N t is the num- 
ber of integration steps per period. 

The integration is carried out computing the coeffi- 
cients c„(ifc + Ai) (tk = kAt) by successive multiplication 
of the vector c n (tk) by a unitary matrix 



c(t k + At) = UQVQ- X c(t k ) 
where U and V are diagonal matrices with 



U n>n = exp 



-ih- 



At 



n 2 - 2n 



A(t) , A{tf 



h 2 



(B2) 



(B3) 



and V n>n = exp(—iAtv n ). 

v n are the eigenvalues of the matrix representation 
M of the operator (1/ft) cos(x). Q is an orthogonal 
(unitary) matrix that transforms M into diagonal form 
(M = QVQ^ 1 ). V is a diagonal matrix with elements 
v n . 

The symmetric matrix M has the structure M n>m = 
(Sm,n+i J rS mt „-i)/(2h). This method permits to decrease 
the number of steps per period [351 ]. Nevertheless, it in- 
volves several matrix products and therefore its efficiency 
is degraded. 

To check our results we use a second integration 
method which is more efficient [29]. We expand again 
over plane waves but this time we write down the ex- 
pansion coefficients as a product of two time dependent 
functions, viz. 



N 



(B4) 



n=—N 



Assuming that <f> n is a solution of the problem 

A(ty" 



h 2 

ihd t (f>„(t) = — 



(B5) 





' 2n f 









±j\( t ') 2 dt^y (B6) 



Then inserting Eq. (|B4j) into Eq. (lAlj) we obtain the set of 
equations for the amplitudes a n (t) 



. 1 

''2h 



^+1 



(B7) 



To obtain the propagator U(T,0), we start with the 
initial states = 5 n .k, where — N < n,k < N and 
integrate over one time period T. Here the column of 
the matrix for the propagator is given by a n {T)<f> n (T) . 
For the computation we have neglected the contribution 
originating from A 2 (t) and the constant term 1/h in the 
potential since they only yield global phase factors. 



APPENDIX C: SYMMETRIES OF THE 
FLOQUET OPERATOR 

Let us first analyze the time reversal symmetry: 

Si) . t ► x ^ x. 

In the momentum representation basis, \n), this trans- 
formation implies several operations including complex 
conjugation, and inversions n — > — n. 

In that case the Floquet operator possesses the prop- 
erty 

U(T,t ) = U(T,t f, 
whose matrix elements obey 



The latter relation holds whenever the transformation is 
done around the symmetric points, namely to = t s . 
Let us analyze the symmetry 

S a :t^t + T/2, x -> -x. 
This corresponds to the transformation 
t^t + T/2, n->-n. 
For the Floquet operator we get 

U(T,0) ee U(T,T/2)U(T/2,0) 



/7(T/2,0)* J 'l7(T/2,0). (CI) 



This implies that 



or, equivalently, 



U ->U 
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APPENDIX D: ASYMPTOTIC CURRENT 



The asymptotic current can be defined as follows 



,/= lim — y (ip(mT)\p\ip(mT))T — 



lim -X>("tfO»T. pi) 
The wave function can be expanded over Floquet states 

as 

|V(*))=£C a |<^>- (D2) 

a 

Similarly, one can expand the Floquet states over the 
plane wave function 

!*«(*)> = £ ). (D3) 
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